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ABSTRA.C!r 


An  ansljrtical  auid  experimental  investigation  of  the  velocity  and 
pressure  distribution  in  the  wake  of  a  body  of  revolution  has  been 
conducted  and  an  equation  was  derived  which  presents  the  velocity  and 
pressure  distribution  of  the  turbulent  wake  as  function  of  sin  experimental 
parameter  ^  .  The  experimental  parameter-^  was  extracted  from  newly 
performed  measurements. 

Satisfactory  agreement  between  theoretlft^l  prediction  and  experi¬ 
mental  results  was  obtained  over  a  range  of  between  four  and  twelve  dia¬ 
meters  downstream  of  the  body. 
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SECTION  1 


INTHODUCTION 

The  velocity  and  pressure  dlstrihution  in  the  turbulent  wake  of 
two  -  and  three-dimensional  bodies  has  been  treated  in  a  number  of  in¬ 
vestigations.  The  subject  is  very  complicated  and  a  perfect  solution 
does  not  seem  to  exist  (Ref.  5). 

For  problems  of  aerodynamic  retardation,  in  which  a  drag-producing 
secondary  body  is  deployed  in  the  wake  of  a  fast  moving  primary  body  in 
order  to  reduce  through  aerodynamic  drag  the  speed  of  both  bodies,  the 
knowledge  of  the  velocity  and  pressure  distribution  in  the  wake  of  the 
primary  body  becomes  a  very  important  problem.  The  final  objective  of 
problems  of  this  nature  is  then  the  determination  of  the  drag  of  retar¬ 
dation  devices  in  undisturbed  flow  and  in  the  wake  of  a  primary  body,  and 
the  drag  of  the  system  consisting  of  primary  and  secondary  body. 

The  analytical  and  experimental  determination  of  the  velocity  and 
pressure  distribution  in  the  turbulent  wake  was  selected  to  be  the  first 
phase  in  this  field  of  investigation,  and  the  approach  and  the  results 
of  these  studies  are  reported  in  the  following. 


Manuscript  released  by  the  authors  on  December  1959  for  publication. 


WADD  TE  60-257 


-1- 


a  turtulent  flow  with  a  varying  velocity  u  which  is  smaller  than  the 
velocity  of  the  moving  body  (Pig.  l).  In  the  other  case  the  body  is  at 
rest  while  the  vindisturbed  fluid  moves  with  the  velocity  Uq,  with  a  vary¬ 
ing  velocity  u'  in  the  wake  of  the  body  (Pig.  2).  Experience  shows  that 
with  increasing  distance  from  the  body  the  area  of  the  disturbed  velo¬ 
city  increases  while  the  difference  between  the  free-stream  and  the  dis¬ 
turbed  velocity  decreases. 


2.2  The  Simplified  Eotmtion  of  Motion 

Let  us  assume  that  a  body  of  revolution  is  immersed  in  a  fluid 
with  the  velocity  Uq  while  the  body  is  held  at  rest  and  that  a  fixed 
longitudinal  axis  is  arranged  which  coincides  with  the  direction  of  the 
undisturbed  flow  (Pig.  2). 

At  a  certain  distance  x  behind  the  body  the  components  of  the  dis¬ 
turbed  velocity  in  the  x-  and  y-directions  are  u'=  u^  -  u  and  v  respec¬ 
tively.  As  is  usual  in  considerations  concerning  boTindary  layer 
phenomena,  one  may  assume  that  the  variation  of  the  velocity  perpendic¬ 
ular  to  the  x-axls  is  large  compared  to  that  along  the-  z-axls,  and  that 
the  pressure  gradient  parallel  to  the  x-axls  may  be  neglected.  Under 
these  circumstances  one  may  establish  with  Swain  (Bef .  3)  l^e  equation  of 
continuity  and  the  approximate  differential  aquation  of  motion  as 
follows: 


^  (yu)  +  {  yv)  0 

y  X  3  y 


(1) 


and 


(u.-u)|M 


f  V 


^  - 
9y 


ay 


'A(— ) 


(2) 
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The  term 


in  equation  (2)  is  the  so-called  "apparent  stress"  or"Eeynolds  stress," 
while  represents  the  "Misch'uneiswec''  (Ref.  1)  Usin^;  the  condition  of 
continuity  and  the  fact  that  the  aerodynamic  drag  of  the  "body  equals  the 
loss  of  momentum  in  the  turbulent  waJce,  one  may  write 


Cp  S  Uo^'=  Zirp  J  (u^+  Uq  u)  y  dy, 


(3) 


in  which  the  term 


^0 


is  the  conventional  presentation,  of  the  aerodynamic  drag  of  the  body. 
Since  the  velocities  u  and  v  are  small  compared  to  u^,  equations  (2)  and 
(3)  can  be  written  as  follows: 

,,  ^  •‘i  _  1  d 

'^o  Tx-  y  dy 

and 

P  2 

^  Cd  s  Uo  =  2irp  /  Uq  u  y  dy.  (5) 

'^0 

In  cylindrical  coordinates  equation  (i)  is  written  as 
p  9 

^  s  Ug*^  =  2  7r^  /  '^o  ^  (5a) 

Replacing  the  term  for  the  apparent  stress  in  the  conventional 
manner,  one  obtains  (Ref.  5) 

^  '  (^)  ~  (%ar  ■  '^in^ 

where  M-  is  an  empirical  constant  and  b  denotes  the  width  of  the  mixing 
zone. 
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With  -•  0  and  the  derivation  completed,  equation  (4)  may  he 


written  in  cylindrical  coordinates  r  and  x  as 


u 


o 


%ax 


/I  ^ 

\r  dr 


172)  • 


(7) 


Assuming  that  at  a  great  distance  from  the  body  the  profiles  of  u 
in  different  cross-sections  of  the  turbulent  waJce  are  similar,  and  using 
the  method  of  analysis  of  dimensions  and  power-laws  (Befs.  3  and  5),  one 
obtains 


r-  7;  1^/3, 
b  =  K  x^^^, 


and  the  velocity  at  the  centerline 


In  these  equations  7]>,  K,  and  A  are  coefficients  of  proportionality. 
The  solution  of  equation  (?)  may  be  assumed  to  be 


u 


A  Uo  3772 


(8) 


in  which  s  needs  to  be  determined. 

The  value  of  the  individual  terms  of  equation  (?)  with  the  assump¬ 
tion  under  (8)  amounts  to 


du 

dx 


^  _  2  A  Up 
dr  X 


s7f 


and 


Sr2  4/3 


(1+2  b7)  ) 
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Substituting  these  expressions  into  equation  (?)  results  in 

(1+817^)  (1  6  K  A  s)  =  0. 

which  provides 


and 


"  772 


82  =  - 


6  K  A 


Equation  (8)  may  now  be  written 


A  uo  /  6  VlK  A 

u  =  — I  C,  e 

,2/3  ' 


-1 


+  C2  e 


The  boundary  conditions  of  the  wake  are 

u  =  0  for  r  =  00  . 
=  '^sax  for  r  =  0 


from  which  the  constants  can  be  evaluated  as 


C,  =  1. 


C2  =  0. 


Therefore,  the  disturbing  velocity  amounts  to 

772 


u  =  u. 


,2/3 


A  6  KK  A 

e 


(9) 


2.3  The  Coefficients  of  A  and  K 

The  velocity  relationship  presented  in  equation  (9)  contains  two  un¬ 
known  coefficients  A  and  K.  For  their  determination  one  has  so  far  merely 
the  drag  momentum  equation  The  second  condition  will  be  taken  from 

experimental  data  given  in  Hef .  10  and  those  shown  in  Fig.  k.  Furthermore, 
the  analytical  treatment  by  Swain  (Ref.  3)  can  be  used  to  derive  another 
condition  for  the  determination  of  A  and  K  for  the  purpose  of  ccmparing 
the  validity  of  the  entire  treatment. 
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For  the  utilization  of  the  experimental  data,  equation  (9)  will  he  trans¬ 


ferred  into  the  form 


u 

%iax 


K  r'^ 
6  VIA 


=  e 


(10) 


with  the  notation 


r'  = 


xl/3 


r 

h 


Experiments  show  that  the  ratio  of  u/u^^^^  varies  very  little  in  the 
vicinity  of  the  centerline,  which  qualifies  this  region  particularly  for 
numerical  evaluation.  Therefore  the  velocity  ratio  at  a  point  with  r'  =  0.1 
shall  he  chosen,  which  amounts  in  accordance  to  Eef.  10  and  Fig.  if  to  approx¬ 
imately  u/u  “0.96.  Therefore  one  may  set 
““  0.01  K 


0.96  *  e 


6  V?-  A 


which  provides 


K  =  24.^  A 


(11) 


Expressing  the  drag-momentvuB  equation  (5a)  in  terms  of  the  velocity 


as  presented  in  (10),  one  obtains 


OO 


K  (  r')^ 

6  A 

r'  dr' 


CpS 

4  TT 


3  A^  K  K 


Substituting  the  previously  obtained  value  for  K  into  the  last  equation 


gives 


CpS  =  12  rr  24.i<8  A^ 


(12) 


or 


A  =  0.103 
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From  equation  (11)  and  (12)  one  finds 


Z  -  2.52  (  CdS  >^  ) 


1/3 


(13) 


Swain  (Ref.  3)  derived  for  the  velocity  distribution  the 
equation 

-2-  .  if  -  if 

'%ax 


(14) 


in  which 


^  =  JL,  ^  JL. 

^  r«  h 


Tq  being  the  value  of  r  at  the  edge  of  the  wake. 

One  may  now  consider  the  velocity  ratio  of  the  particular  point  at 
which  u  =  Umax  call  the  related  coordinate 


i  -c 


This  coordinate  can  numerically  be  determined  from  the  relationship 


'^mai  3/2  ,^2 

u  =  __  =  W  m  - 


which  provides 


f m  =  0-293 


2/3 


Equation  (lO)  can  be  transformed  by  means  of  rj^  as  the  coordinate  at  which 
%ax 


u  = 


in  the  same  manner  and  one  obtains 


u 


^ax 

2 


Z  rm^ 

'  6 

Umax  ® 
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From  this  one  finds 


Pni  =  (  4.158 


ye  A  .1/2 
K 


If  one  postulates  that  both  systems  are  correct,  the  Talues  ^  211 
and  rjQ  must  he  eqital.  Consequently, 


(  4.158  =  (  0.293 


from  which  follows 


A  =  0.0468  . 


Suhstitutin^  A  into  (12)  proTides 


S  =  12  TT  (  0.0486)^  -|- 


whence 


and 


A  -  2.3  (  Cp  S 
A  =  0.107  {  Cp  S 


(17) 


The  values  of  A  and  A  which  were  obtained  from  experimental  data, 
shown  under  (12)  and  (I3),  a^nd  the  related  data  derived  from  Swain's 
equation,  presented  above,  agree  very  well.  In  view  of  the  fact  the  results 
from  the  experimental  data  as  well  as  the  ones  derived  from  Swain's 
equation  include  certain  assumptions,  it  appears  to  be  Justified  to 
establish  average  values  for  A  and  A  based  on  the  two  independent  methods. 
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In  this  manner  one  can  establlah 


A  =  0.105  (  Cjj  S  /«? 
K  =  2.41  (  Cjj  S 


These  values  shall  he  used  In  the  further  treatment  of  the  woke  problem. 

The  velocity  equation  (9)  presented  with  7^  -  ■  y-  and  with  A  and  Z 

X  ^ 

from  above  may  now  be  written  in  the  following  form: 


u _  _  O.105  / 

'^o  x2/3 


S  1'^'  ) 


1.525  (  Cj3  S>tx)2/3 


(18) 


SECTION  3 

NUMEKICAL  COMPAfilSON 

The  equations  (ih)  amd  (18)  represent  the  local  velocity  In  the  tur- 
b\ilent  wake.  Judging  from  their  structure,  It  appears  that  equation  (18) 
might  offer  an  easier  numerical  process  and  it  appears  to  be  Justified  to 
pursue  Its  exploitation.  The  first  point  In  this  study  shall  be  a  com¬ 
parison  of  the  predictions  of  available  theoretical  methods  with  experi¬ 
mental  resxilts.  For  this  purpose,  one  may  compare  the  predictions  re¬ 
presented  by  equation  (18)  with  those  given  by  Swain's  equation  (14) , 
Goldstein's  results  obtained  by  meane  of  "modified  vortlcity  theory" 

(Bef.  9),  £uid  the  mean  experimental  results  obtained  by  Hall  and  Hlslop 
(Eef.  10).  For  a  convenient  comparison  one  may  Introduce  the  quantity 
or  Cm  defined  before,  and  derive  from  equations  (l4)  and  (18) 
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u 


2 


(19) 


0.293 


and 


(20) 


respectively.  Equation*  (19)  and  (20)  represent  now  the  velocity  ratio 


as  functions  of  =  z.  The  same  process  can  also  he  applied  to 

m  ^  ID 

the  findings  given  hy  Goldstein's  (Hef.  9)  a&d  the  experimental  results  by 
Hall  and  Hislop  (Hef.  lO).  In  the  manner  the  results  of  Table  1  were 
calculated.  Figure  3  illustrates  the  same  results  and  it  c£ui  be  seen  that 
equation  (20)  agrees  with  the  experimental  values  best  in  the  central 
portion,  while  equation  (19)  appears  to  be  a  better  approximation  in  the 
outer  portion  of  the  wake. 


3.1  Eevised  Equations  for  the  Velocity  Distribution  behind  a  BocLv  of 
Revolution 

For  practical  cases  the  drag  coefficient  is  usually  known  and  the 

ITT? 

related  area  is  generally  given  by  S  =  ^  ,  where  Die  the  diameter  of 

the  body  of  revolution.  In  such  a  case,  it  is  convenient  to  present 
equation  (18)  in  terms  of  Cp  and  D  and  one  obtains 
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1/3 


(21) 


u  ^  0.107 

^ \  ^  / 


_  0.415  (r*)^ 

ji|/3  (CD7r>t)2/3 

e 


where  r*  = 


Thus,  the  Telocity  In  the  wake  of  the  hody  of  revolution  which  moves 
with  a  velocity  through  the  fluid  is  given  hy  equation  (21).  If  the 
fluid  moves  with  the  velocity  u^  past  the  fixed  hody,  the  velocity  of  the 
fluid  in  the  wake  u'  equals 


u'  -  u^  -  u  , 

and  may  he  expressed  hy  means  of  equation  (21)  as  follows: 


0.415  (r*)' 


(CdTTK) 


2/3 
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TABLE  I—  ANALYTICAL  AND  EXPERIMENTAL 
VELOCITY  DISTRIBUTION 
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AND  EXPERIMENTAL  STUDIES 


SECTION  4 


THE  EXPERIMEKTAI.  PABAMETER 

The  equations  (21)  and  (22)  still  contain  the  empirical  parameter  vt. 
which  for  any  practical  solution  needs  to  he  known.  For  the  purpose  of 
this  study  and  for  a  more  experimental  approach  to  the  principal  problem 
of  aerodynamic  retardation,  the  pressure  distribution  in  the  wake  of  a 
cylindrical  body  with  an  ogive  nose  and  blunt  tall  was  experimentally 
determined.  Some  results  of  this  study  are  presented  in  Fig.  4  and 
Table  2,  which  are  abstracts  of  Ref.  12.  An  attempt  will  be  made  to 
correlate  the  exiierlments  with  the  analytical  prediction  of  equation  (21) 
and  various  values  of  the  parameter  Vt  , 

4.1  Experimental  Pressure  Dietribution 

The  pressure  distribution  of  the  wake  shall  be  presented  by  means  of 
the  pressure  coefficient,  which  is  defined  as  the  difference  between  the 
total  and  the  static  pressure  in  the  wake  divided  by  the  dynamic  pressure 
of  the  free  stresun.  For  incompressible  flow  this  fraction  is  Identical 
to  the  ratio  of  the  square  of  the  velocities  in  the  wake  and  in  the  free 
stream.  In  view  of  the  derivation  above  one  may  set 


which  can  be  expressed  in  terns  of  equation  (22)  and  one  obtains 

-,2. 

0-415  (r*)^ _ ^  1 

0.107  / 

X. 2/3^4  ^2 

Di 


.1/3 


/ir/3 


(Cd  tt  «.) 


2/3 


(23) 
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TABLE  2—  RESULTS  OF  EXPERIMENTAL  DETERMINATION  OF  THE 
PRESSURE  DISTRIBUTION  IN  THE  WAKE  OF  A  BODY  OF 
REVOLUTION  (M~0.2,  =2.74  ^  lO"",  Cp=0.35) 


li. 
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PRESSURE  DISTRIBUTION  I 
BODY  OF  REVOLUTION 


In  the  equations  (21),  (22),  and  (23)  the  velocity  ratio  ^ —  depends. 

for  given  Cj,  ,  ^  ,  and  r*  values,  only  upon  the  parameter  X.  ,  which  can  he 

determined  from  any  point  at  which  —  is  known.  From  Fig.  4,  it  is  evident 

'lo 

that  M.  will  vary  depending  on  the  location  of  the  control  section.  Other¬ 
wise  it  is  convenient  to  choose  points  at  the  longitudinal  sucis  as  refer¬ 
ence  points.  With  r*  =  0  and  =  O.35,  equation  (2l)  becomes 


u  ^  Q.1C7  /  . 


from  which  H  is  found  to  he 


0.0176 
7  /u\3/2 


With  equation  (24)  and  the  experimental,  data  given  in  Fig.  4  one  cam 
determine  for  each  2  value  a  related  value  of  hL  as  shown  in  Table  3 
amd  Fig.  5. 


TABLE  3 


0.525 

0.725 

0.775 

0.800 

O.8I5 

0.630 

0.275 

0.149 

0.120 

0.106 

0.097 

0.089 

0.00372 

0.00584 

0.00497 

0.00406 

0.00339 

0.00305 

0.0610 

0.0764 

0.0705 

0 . 0633 

0.0582 

0.0552 
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X  VS.  X/D  AS  DETERMINED  FROM  EXPERIMENTAL  DATA 


The  average  value  of  ^  from  Table  3  amounts  to  K  -0.0633,  which 

corresponds  to  Cp  at  ?  =.  8.  Substituting  ^  -0.0633  and  =  0.35 
•^0  D 

into  equation  (2l)  provides 


Corresponding  to  this  relationship,  the  pressure  distribution  was  calcu¬ 
lated  and  the  results  are  shown  in  Table  4  and  Pig.  6. 

An  inspection  indicates  good  agreement  between  experimental  and 

analytical  data  at  distances  of  -  =  6. 

D 

In  an  attempt  to  obtain  a  better  agreement  between  analytical 
prediction  and  experimental  results  at  shorter  distances  from  the  body, 
the  velocity  and  pressure  distribution  was  calculated  from  equations  (22) 
and  (23)  using  it  -•  0.0764  corresponding  to  an  ^  value  of  4  (  See 
Fig.  5  or  Table  3). 

The  numerical  forms  of  equations  (2l)  and  (23)  amount  to 


Table  5  and  Fig.  7  show  the  results  of  these  calculations.  One  observes 
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FIG.  6  EXPERIMENTAL  AND  ANALYTICAL  DATA  FOR  A  BODY  OF  REVOLU¬ 
TION  (BASED  ON  ^=0.0633  RELATED  TO  Cp  AT  X/D=  8) 


TABLE  5 —  RESULTS  OF  ANALYTICAL  DETERMINATION  OF  THE 
PRESSURE  DISTRIBUTION  IN  THE  WAKE  OF  A  BODY  OF  REVOLUTION 
(  BASE->D  ON  9^=  0,0764  RELATED  TO  Cp  AT  X/O  =  4  ) 
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FIG. 7  EXPERIMENTAL  AND  ANALYTICAL  DATA  FOR  A  BODY  OF  REVOLUTION 
(BASED  ON  ^=0.0764  RELATED  TO  Cp  AT  X/D=4) 


now  an  Improved  agreement  Isetween  theory  and  experiment  in  the  central 
portion  of  the  wake  In  the  range  of  x/D  =  4  to  x/D  =  8,  whereas  the 
results  In  the  more  outward  portion  differ  as  before.  Also,  the 
results  for  x/D  >8  do  not  agree  as  well  as  before. 

The  above  procedure  was  again  repeated  using  the  value  >€.  =0.0552, 
related  to  x/D  =  12.  Table  6  and  Pig  8  show  the  results  of  these  calcula¬ 
tions.  One  obseives  that  the  analytical  and  experimental  results  agree 
well  only  for  i/D  =  10. 


SECTION  5 

CONCLUSION. 

An  equation  was  derived  which  presents  the  velocity  and  pressure 
distribution  of  the  turbulent  wake  as  function  of  an  empirical  parameter 

A  study  revealed  that  the  predictions  resulting  from  the  newly 
derived  formula  agree  satisfactorily  with  those  from  known  theories  and 
experimental  data. 

The  empirical  parameter  ^  was  extracted  from  newly  performed  experi¬ 
ments  and  It  was  found  that  depends  strongly  upon  the  location  of  the 
control  section  In  relation  to  the  rear  end  of  the  wake  producing  body. 

Palrly  good  agreement  between  theoretical  prediction  and  experimental 
results  was  obtained  by  means  of  am  average  value  of  =  0.0633  over  a 
range  of  J 

Por  better  agreement  at  a  particular  location,  especially  closer  to 
the  primary  body  or  at  great  distance  from  It,  specific  ^  -values  are 
recommended  which  are  presented  as  function  of  the  dimensionless  distance 

The  numerical  values  of  ^  presented  in  this  study  are  related  to 
a  wake  producing  body  having  drag  coefficients  of  =  0.35*  It  may  be 
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t)  IH 


expected  that  the  X.  -value*  at  function  of  x/l»  depend  slightly  on  C  of 
the  primary  hody. 
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0.737  I  0.780  0.309  O.83O 
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TABLE  6—  RESULTS  OF  ANALYTICAL  DETERMINATION  OF  THE  PRESSURE 
DISTRIBUTION  IN  THE  WAKE  OF  A  BODY  OF  REVOLUTION 
(BASED  ON  3^=0.0552  RELATED  TO  Cp  AT  X/D=I2) 
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